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A RADIOSONDE THERMAL SENSOR TECHNIQUE FOR 
MEASUREMENT OF ATMOSPHERIC TURBULENCE 


Jack L. Bufton 
Goddard Space. Flight Center 


INTRODUCTION 

In the past few years increased interest in propagation effects in the free atmosphere at mi- 
crowave and optical frequencies has encouraged the study ot microscale turbulence. Even 
in clear air the density inhomogeneities resulting from turbulent mixing tend to scatter in- 
cident radiation. The resulting phase and amplitude degradation of the propagated radiation 
can lead to serious problems for radar and communication links. Theory relating turbulence 
structure to propagation statistics is generally available for the interesting propagation paths 
but, until recently, little turbulence data had been gathered. This lack of data is especially 
severe for the atmosphere beyond a few tens of meters above ground level. Early attempts 
to collect these data concentrated on profiles and the spectra of temperature, wind velocity, 
energy dissipation rate, etc. Turbulence theory was then required to convert the data into 
a strength factor for prediction of propagation statistics. Recently, remote sensing techniques 
of radar backscatter and optical correlation have been used or proposed to infer the necessary 
profiles of the turbulence strength factor. Other recent efforts have concentrated on direct 
point measurement of this strength by using temperature structure as a tracer for turbulence 
structure. Typical results of these efforts are the data reported by Tsvang 1 and Lawrence, 
Ochs, and Clifford 2 on the strength of microthermal fluctuations in the free atmosphere to 
heights of a few kilometers. The work reported here combines the thermal sensor technology 
for microthermal measurements with radiosonde balloon systems. The result is an extension 
of turbulence sensing for the first time to heights up to 25 km above sea level. The result is 
a complete vertical profile of microscale turbulence structure for use in analysis of space-to- 
earth and earth-to-space propagation links. The measurement system is here referred to as 
the thermosonde (thermal sensor radiosonde) system. Thermosonde data provide the quanti- 
tative link between theory and propagation statistics. The data may also be used to infer the 
large scale turbulence structure known as CAT (clear air turbulence), which adversely affects 
aircraft operations. 

In the thermosonde system a pair of high resolution thermal sensors and the associated 
electronics report temperature statistics by telemetry during both balloon ascent and descent. 
The prototype payload instrumentation was tested during October 1970 and again during 
September 1971, in conjunction with an earth-to-balloon laser propagation experiment at 
White Sands Missile Range, New Mexico. The thermosonde payload was reeled down from 



the main balloon payload to gather turbulence data for analysis of the laser propagation re- 
suits. Details of the prototype payload plus examples of the turbulence profiles are reported 
by Minott, Bufton,and Fitzmaurice 3 and by Bufton et al. 4 The production thermosonde 
payload is an improved light-weight, low-cost version of the prototype. Because of these 
improvements it is suitable for use with standard radiosonde balloon systems. Both the pro- 
totype and the final payload version were developed and produced by GTE-Sylvania under 
contract to the Laser Technology Branch, GSFC. 

This report will present specifications and operation details of the final thermosonde payload. 
The turbulence profile data collected will be presented and discussed in light of the tempera- 
ture and wind velocity fields for the area of each flight. Applications of the turbulence pro- 
files to analysis of optical propagation results will also be discussed in detail. In particular 
the results of near simultaneous stellar observations and thermosonde flights are compared. 
The major motivation for this comparison as for the thermosonde development program it- 
self was the desire to monitor atmospheric turbulence effects at high altitude on laser com- 
munication links. Progress toward this goal and other applications of the thermosonde 
data are discussed in the comments section. 


DESCRIPTION OF THERMOSONDE SYSTEM 
Theory of Measurement 

The fluctuations in refractive index that result from turbulent mixing are quite random in 
nature and are best described by statistics. The fundamental statistical quantity is the structure 
function and its strength factor, the structure coefficient. 

Turbulence strength and distribution enters optical propagation theory through C^(h), the 
refractive-index structure coefficient. Propagation statistics invariably contain integrals over 
this turbulence strength factor as it varies with altitude h. For optical frequencies, (h) 
depends almost entirely on microscale temperature fluctuations specified by the temperature 
structure coefficient, C 2 (h). The connection at a wavelength of 0.5 pm is provided by 



80 P(h) 
T(h) 2 



( 1 ) 


where P(h) = pressure (mb) and T(h) = temperature (K). The thermosonde system pro- 
vides C^(h) data and thus Cj* (h) data by measuring mean square temperature fluctuations as 
a function of altitude and using the widely accepted Obukhov^ -Kolmogorov^ turbulence 
theory, 
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This theory requires isotropic turbulence.structure, and is applicable for turbulence element 
sizes (r) larger than 1^ the inner scale, and smaller than l Q ,the outer scale. Turbulence elements 
between these scale sizes comprise the inertial subrange. Energy input to the turbulence 
process through structures of dimension 1 Q cascades down adiabatically through the inertial 
subrange until it is lost as heat in structure of dimension less than 1.. Tsvang 1 and others have 
indicated that the existence of an inertial subrange and the validity of the 2/3-power law 
are well supported by field measurement. The assumption of isotropy indicated above must be 
modified to one of local isotropy, since turbulence strength and scale sizes 1 Q and l t are expected 
to vary with altitude. Thus in equation (2),altitude, h, is indicated as a parameter of the 
measurement. The scale size 1 4 is on the order of a few cm near the surface and increases to 
tens of cm in the stratosphere. The outer scale 1 is thought to be roughly equal to altitude 
above the surface. The scale sizes which most affect the optical propagation generally fall 
within these inertial subrange limits at all altitudes and are weighted most heavily towards 
those near 1 j The end result of the microthermal measurements when expressed in terms 
of refractive index statistics is . Its variations with altitude constitute a vertical profile of 
the optical strength of turbulence. 

Details and Specifications of the Thermosonde System 

The thermosonde payload is used to measure root mean square (rms) temperature difference 
between two fine platinum wires, and to send the result over telemetry to a ground receiver 
station. Each platinum wire probe is 2.5 pm in diameter and about 6 mm long. They have a 
linear resistance-temperature coefficient, and respond to an increase in atmospheric tempera- 
ture with an increase in resistance. For a temperature change, AT, of 1.0 K a typical re- 
sistance change, AR, is near 0.3 ohm. Resistance changes accurately follow temperature 
changes at frequencies from dc to near 1 kHz. The two probes are mounted parallel to each 
other and separated by 1 .0 m in the horizontal plane. Each fine wire probe is soldered to two 
stiff wires about 8 cm long, which are in turn soldered to a printed circuit card. When the 
two cards are attached to polyurethane support arms the distance between probes is 1 m. 

The probes are connected as two arms of a Wheatstone bridge. This permits a sensitive meas- 
ure of temperature difference. The bridge is driven by a 3-kHz oscillator. The bridge output 
carrier signal is amplitude-modulated by the temperature difference signal. Figure 1 
is a block diagram of the bridge temperature measurement scheme, the electronics for com- 
puting mean square temperature difference, and the telemetry transmitter. 

The bridge output voltage is amplified, bandpass filtered and synchronously detected. After 
demodulation, the signal is low-pass filtered to 1 kHz, then sent through an rms operator mod- 
ule with a 5-second RC output time constant. The result is a voltage proportional to root mean 
square temperature difference. Calibration of this output is achieved before flight by tem- 
porarily inserting a reference resistor into the bridge network. A feedback system senses the 
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Figure 1. Thermosonde Signal Processing Electronics 


4 






average AT voltage and selects unity gain or 13X gain at the input to the nils module. Gain 
switching effectively increases the dynamic range of rms measurement to 44 dB from 0.004 K 
(system noise level) to 0.676 K. The system noise level yields a (AT) 2 of 1.6 X I O' 5 K 2 and 
from equation (2) with r = lm, Cj = 1 .6 X 10' 5 C 2 m‘ 2/3 . This in turn with equation (1 ) 
implies a typical noise level forC^(h) of 1.6 X 10‘ 17 m ' 2 ^ 3 at surface level; 2.9 X 10' 18 m" 2 ^ 
at h = 10 km and 1.4 X 10' 19 m‘ 2/3 at h = 20 km. Another feedback loop corrects for bridge 
imbalance over an RC time constant of 40 seconds. Long term imbalance could result from 
small differences between resistance-temperature coefficients of the probes, and the typical 
ambient temperature changes from 290 K near the surface to 200 K at 10 km altitude. Be- 
cause of this feedback loop the low end response to temperature changes does not extend 
entirely to dc. The net accuracy for rms computation is within ±10 percent for the extremes 
of turbulence and ambient temperature indicated above. 

Rms output is sent along with a dc voltage indicative of gain scale to the telemetry trans- 
mitter. The telemetry system is a standard Weather Service frequency-modulated 1680 MHz 
system with transmitter modified to have a linear response to input voltage over a 46-dB 
dynamic range and to delete the standard temperature and relative humidity data. Rms 
temperature fluctuation data arc sent back continuously as the payload ascends and descends 
except for brief interruptions controlled by a pressure activated commutator during which 
scale information is transmitted. 

The entire telemetry system was calibrated by applying a voltage equal to AT max of 0.676 K 
and one equal to AT min of 0.004 K to the transmitter. At the receiver, the output, which 
appears on chart paper, was adjusted to properly record these high and low limits. The re- 
ceived rms temperature difference data are squared to obtain C 2 (h) and thus (h). 

Figure 2 is a photograph of an actual flight-configured thermosonde payload. The two 
thermal sensor probe supports, visible at either end of the horizontal arm, are separated by 
one meter. The rms signal processing electronics are located on a printed circuit board, 
which is approximately 20 cm square. The board is mounted vertically in the central portion 
of the payload. Pieces of polyurethane thermal insulation 3.8 cm thick enclose and protect 
the electronics, and form the probe support arm. Because of the insulation, temperature 
extremes inside the payload are confined to 273 ± 25 K for a typical flight, A typical probe 
support is composed of two 16-gauge bus wires, soldered to a printed circuit card strip. Hori- 
zontal probe mounting plus the thin rectangular shape of the payload tend to align the 
probes perpendicular to any air flow. The radiosonde modulator/ transmitter is housed in 
the rear section of the payload. A battery supplying the + 1 1 5V and + 6V needed to operate 
this part of the system is also located in the rear section. The battery supplying ± 1 8V to 
the rms electronics and probe circuits is strapped to the top of the payload. Both batteries 
are activated by immersion in water just prior to launch, and provide approximately 3 hours 
of continuous operation. 
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Figure 2. Thermosonde Flight Payload 



'Hie thermosonde payload is mounted approximately 5 meters below the parachute and 
approximately 55 meters below the radiosonde balloon. This is done to minimize balloon 
wake effects on the turbulence data during the ascent phase. Total payload weight is just 
under 2600 grams. Of that total, the major parts are activated batteries (920 gm), processing 
electronics (454 gm), and radiosonde transmitter (480 gm). Table 1 is a summary of the 
thermosonde performance specifications. The thermosondc system final report prepared 
by Titterton, Mallery, and Arken 7 at GTE-Sylvania contains a detailed collection of pay- 
load specifications. 

Table 1 

Thermosonde System Specifications 

Parameter measured: <^(AT) 2 ^> / ^ Cj. 

Measurement time constant (r): 5 seconds 

Measurement range: 0.004 K ^<^(AT) 2 ^> ^ 0.676 K 

J 

AT noise level: 0.004 K maximum 

AT saturation level: 2.0 K 

Frequency range of AT measurement: 0.16 Hz to 1.0 kHz 

diameter: 2.5 pm 
length: 6 mm 

Thermal sensor probe material: platinum 

« separation: 1 m 

Balloon ascent rate: 5ms' 1 

Payload descent rate: 10 ms" 1 

Maximum altitude: 25 km 

Turbulence structure vertical resolution: i m (ascent) 1 

» 50 m (descent ) ) 

Payload weight: 2600 gm 

Telemetry frequency: 1 680 MHz 

Battery power requirements: 1 ^ V at 175 mA +£ v at 300 mA 

+ 115 V at 45 mA 

Typical mission duration: 2.5 hours 
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A typical thermosondc mission profile would include the following: 

1. Battery activation and telemetry calibration before launch. 

2. Ascent at 5ms' 1 to an altitude above 20 km. 

3. Separation of payload and balloon by firing of explosive cutter at a predetermined 
time interval after launch. 

4. Payload descent by parachute at 10 m s" 1 until surface impact. 

During the flight, the balloon and payload follow a trajectory dictated by the mean wind 
structure. M icrothermal turbulence data are sent by telemetry on both ascent and descent. 

The product of balloon vertical speed and the rms time constant (5 seconds) yields a turbulence 
element resolution of 25 m on ascent and 50 m on descent. Thus, the reported rms temper- 
ature data at any instant of time reflects an average over the last 25 or 50 meters. 

Error Mechanisms 

There are several readily identifiable sources of error inherent in the thermosonde turbulence 
measurement technique. As mentioned above the system noise level is at or below AT = 

0.004 K. This is the point where signal-to-noise ratio (SNR) - l .0. Data accuracy increases 
for AT above this value in proportion to signal strength. For AT = 0.4 K, SNR = 40 dB and 
the measurement error is on the order of 1 percent. Additional error bounds of ±10 percent 
were noted above for the rms computation. It is estimated that telemetry errors and reading 
errors at the ground station are considerably less than this ±10 percent. After received data 
is squared to obtain C|(h) and thus, (h), the error is on the order ±20 percent. Since the 
turbulence coefficients C|(h) and C^(h) are expected to vary over several orders of magnitude, 
this error is not of major proportions. 

As noted above, the turbulence clement resolution was limited to 25 m on ascent and 
50 m on descent. There is strong evidence from radar probing that the minimum vertical 
extent of some turbulence layers may be only a few meters. Any such fine structure in the 
turbulence profile would not be apparent in detail in the thermal sensor data, but would of 
course contribute to the 5-second integrated value for Cj.(h). 

The temperature measurement of the sensors may be in error due to a variety of factors in 
the heat exchange of each sensor with the air, radiation field, and support structures. Heat 
exchange errors would be significant only when they have different relative effects on each 
sensor. This is because temperature difference, not absolute temperature, is the parameter 
measured. A likely candidate for differential error is daytime solar radiation. 

During the ascent portion of the flight, air flow over the balloon and payload may break into 
turbulence, and thus cause the thermal sensors to report higher turbulence strengths. The 
most serious wake effects should occur at high altitude, where the temperature differential 
between balloon and atmosphere contributes thermal eddies to the wake. For this reason 
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descent data will be used to predict the optical effects. Unfortunately the descent profiles 
do not include data in the first 1 or 2 km above mean sea level (MSL) because ot telemetry 
dropouts. Thus it will be necessary to use some ascent data to have a complete profile for 
comparison with the optical data. 

Turbulence non-stationarity also affects the quality of the thermal sensor data, and in particu- 
lar its usefulness in computing optical propagation statistics. Although it is helpful that 
major turbulence layers and features often are visible on both ascent and descent profiles, 
even though the profiles are separated by tens of minutes and tens of kilometers, the fact 
remains that thermosonde data result from point measurements of structure changing in 
both time and space. For effective correlations with propagation data it is essential that 
the turbulence data be acquired simultaneously and that the balloon flight path be near 
the propagation line of sight. 

A more serious source of error may be the possible departure of the turbulence spatial spec- 
trum from the classical Obukhov-Kolmogorov (-5/3) power law for wave number. This theory 
implies the existence of an inertial subrange ot turbulence element sizes for which the turbu- 
lence is assumed to be, at least locally, isotropic and homogeneous. For well developed tur- 
bulences, a large amount of evidence supports this description. However, for a developing 
turbulence and for the breaking wave motion characteristic of upper altitude turbulence 
structure, another power law may apply. In addition, the limited vertical extent of some 
observed structure indicates a small outer scale, l 0 , and truncation of the inertial subrange. 

The major effect of this error source would be to alter the quantitative prediction of optical 
effects. 

Despite these rather serious qualifications of the observed data, the tact remains that the 
thermal sensor data are a fairly accurate point measurement of the temperature structure func- 
tion for the spatial separation of one meter. As such, the data serve as at least a relative indica- 
tion of the optical strength of turbulence, and a quite accurate indication of turbulence location. 

MICROTHERMAL TURBULENCE DATA 
Thermosonde Data Collection 

Turbulence data were obtained with the prototype thermosonde payload during October 
1970 and again during September 1971, as mentioned in the introduction. These flights 
were conducted as part of the Balloon Atmospheric Propagation Experiment (BAPE) at White 
Sands Missile Range (WSMR), New Mexico. One low-altitude turbulence profile was obtained 
during BAPE I (1970). This result, plus an illustration of the application of turbulence data 
to the laser propagation problem, is contained in the BAPE T final report. 3 

During BAPE II (1971) four turbulence profiles up to 15 km above MSL were obtained. 

All flights were from the runway of Holloman Air Force Base near Alamogordo, N. M. 
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The base is located on the desert floor, 1.2 km above MSL, and is bordered 50 km to the 
west by the San Andres mountains and 20 km to the east by the Sacramento mountains. 

Both mountain chains run north-south and rise 1 to 2 km above the desert floor. 

The vehicle was a polyethylene balloon of 40-m diameter, which carried aloft a 450-kg optical- 
propagation experiment payload. The complete thermal sensor system weighed less than 
45 kg, and was lowered on a nylon cord 300 m beneath the main payload when the balloon 
reached an altitude of 1 km. The purpose of lowering it was to minimize the balloon-wake 
effects on the thermal sensor measurement of turbulence as the balloon ascended 5 m s' 1 . 

Alt launches were at dawn, in order to have minimum surface winds at the launch site. Thus, 
all data are for turbulent conditions before, during, and up to one hour after sunrise. 

Turbulence data from BAPE II are reported by Bufton et al. 4 Unfortunately, no laser propaga- 
tion data were available for comparison with the thermosonde results. Instead, major emphasis 
was placed on correlation of turbulence activity with indications of turbulence suggested by 
standard meteorological profiles of ambient temperature and wind velocity. These correla- 
tions are reported in detail by Bufton 8 and described in the next section. 

The production version thermosonde payload was first test flown from Chico, California by 
GTE-Sylvania as part of the development program. Thermal sensor data were successfully 
recorded, and the major part of the payload was recovered intact after it had been floating in 
the Pacific Ocean for several days. Subsequent flights of the production model, from which 
data have been recorded and reduced, were launched at the Goddard Optical Research Facility 
(GORF), located near Beltsville, Maryland. This launch site has an elevation of about 50 m 
above mean sea level, with local terrain characteristic of rolling, partially wooded coastal plain. 
The data presented here were obtained under nightime or daytime conditions on several oc- 
casions from May to October, 1972. Most launches were conducted just after dusk or dawn 
in order to minimize winds at launch. About 60 percent of the flight payloads were recovered, 
returned to GSFC for refurbishment, and flown again. 

Both the prototype flights at WSMR and the production version thermosonde flights at GORF 
were accompanied by standard radiosonde measures of upper air temperature and wind velocity 
profiles. At WSMR the radiosonde launches were within a two km horizontal distance and 
within one hour of the main balloon launch. The GORF flights were conducted within 1 to 2 
hours of the National Weather Service radiosonde flight from Sterling, Virginia. The Sterling 
launch site is some 50 km due west of GORF. In both cases the thermosonde and standard 
radiosonde balloons followed similar trajectories dictated by the mean wind structure. Thus, 
both types of sensors sampled nearly the same volume of air. 

Data Reduction 

The baseband output of the 1680-MHz ground-based receiver formed the basic thermosonde 
data set. This output was a voltage proportional to rms temperature fluctuations and was 
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recorded continuously on a single-channel chart recorder. The event marker of this recorder 
was connected to the slow code output of a time code generator. This provided a time base 
(in universal time) which was updated every minute and permanently recorded alongside the 
thermosonde data. The chart record was 1 1.4mm wide with 50 equal divisions. Deflections 
on the chart record were linear with applied voltage,and by virtue of the payload electronics 
and modifications to the radiosonde modulator/transmitter were linear with rms temperature 
fluctuations. Final calibration about two minutes prior to launch was used to establish rms 
temperature reference levels on the chart record. Two references were selected by alternately 
shorting wires on the thermosonde payload. The low reference was used to simulate the 
thermal sensor noise level. For this the chart recorder pen was adjusted to 0 divisions. High 
reference simulated the maximum output AT rms = 0.676 K. Chart recorder gain was ad- 
justed so that high reference occurred at 45 divisions (90 percent of full scale). 

During the thermosonde flight, thermal sensor data were recorded continuously except for 
occasional interruptions by processing electronics gain scale information. These appeared 
as constant chart recorder levels, 25 percent of full scale for low gain and 60 percent of full 
scale for high gain. Scale indications last for several seconds, and were selected every 15 con- 
tacts on the pressure activated commutator of the radiosonde modulator/transmitter. The com- 
mutator calibration curve for contact number versus pressure (nib) was supplied with each unit. 
From this curve and the time history of contact number from the chart record, it was possible 
to plot pressure (mb) of the balloon payload versus time. This was converted into payload 
altitude (MSL) versus time by means of the standard radiosonde record of pressure versus 
altitude. 

A thermosonde data point was selected for each 15-second interval of the chart record by 
visually averaging the raw data values and trends over that interval. Thus, one data point 
was available for each 75-m vertical section of the atmosphere on ascent, and 150-m section 
on descent. The 1 5-second interval included three time constants of the rms operation time 
constant aboard the payload. The selection of average points in effect digitized the raw analog 
chart record. The standard radiosonde data on upper air temperature and wind velocity 
were also digitized, but with coarser resolution (approx. 400 to 500 m). Both data sets were 
entered into a digital computer programmed to convert recorder divisions into rms tempera- 
ture fluctuations, Cj(h) and eventually C^,(h), using the altitude-time curve, the temperature 
and pressure profiles, and the formulae of equations (1) and (2). The computer was also pro- 
grammed to output integrals of the turbulence structure, compute optical propagation 
statistics, and permit rapid and accurate plotting of the input data and program results. Data 
presented below were all plotted in this manner. 

Thermosonde Turbulence Data 

Figures 3 through 9 contain data from all thermosonde flights launched between 27 April 
1972 and 10 October 1972, for which the raw thermal sensor data were reduced to profiles 
of Cj^ (h). Thirteen flights were conducted in this series: Of the six flights for which no 
data are presented, three had in-flight problems with the processing electronics or the 
thermal sensors, two were system checkout flights with no data on descent, and one had a 
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radiosonde telemetry failure. Both ascent and descent (h) profiles for each of the 
remaining seven flights are plotted in the figures. Ascent data for all flights indicated stronger 
turbulence with more layers than the corresponding descent profiles. The consistency of 
this effect suggests the presence of increased turbulence due to balloon wake on ascent. 

During descent by parachute there are no wake effects since the thermal sensors are the 
lead element. 

The major characteristic of both ascent and descent (h) profiles is the irregular pattern 
of multiple spikes. Some of the spikes, or Cj*, (h) maxima, are grouped together to form 
larger structures. Individual spikes, however, have a small vertical extent, often no larger 
than the thermosonde resolution of 25 to 50 m. Despite the ascent data bias due to wake 
effects, there is often good correlation between the larger structures visible on both ascent 
and descent profiles. Good examples are found at 10-1 1 km for flight 2, 9-10 km for flight 
7, and 9-1 1 km for flight 9. The presence of correlation between individual spikes ob- 
served on both ascent and descent is somewhat rarer but still valid. Good examples are 
found at 7.3 km for flight 2, 2.2 km for flight 7, and 9.2 km for flight 9. The presence 
of even limited correlations argues for validity of the thermosonde measurement technique. 

It should also be realized that ascent and descent data sampling are typically separated by 
tens of minutes in time and kilometers in space. Thus, the degree to which correlation is 
present for turbulence structure at the same altitude argues for stationarity of the turbulence 
process and a large horizontal extent for a given structure. The combination of a small 
vertical extent and large horizontal extent suggests that the observed turbulence is in the 
form of layers or horizontal sheets. Lawrence, Ochs, and Clifford 2 noted similar stratifica- 
tion in microthermal measurements performed with aircraft up to an altitude of 3 km above 
ground level. Additional evidence for a multiple layered turbulence structure is contained 
in the radar backscatter measurements performed by Hardy and Katz 9 and Doviak, Goldhirsh, 
and Miller. 10 Their data showed frequent layers at a variety of altitudes up to and including 
the tropopause region. The layers of clear air turbulence seemed to be formed from wavelike 
structures that appear in horizontally stratified sheets. Measured wave amplitudes averaged 
around 500 m, with wavelengths on the order of kilometers. Higher resolution radar studies 
in the lower troposphere (near 300 m altitude) reported by Gossard, Jensen, and Richter 11 
revealed turbulence structure as thin as several meters. Thus, the peaks in C|(z) visible in 
the present data are quite likely the result of a vertical cross section of thin horizontal tur- 
bulence sheets. 

The strongest layers of turbulence of order (^(h) = 10' 14 m' 2 ^ 3 were generally confined to 
the surface boundary layer of about 10 m thickness. Other strong layers of order (h) = 
10- 16 m- 2/3 were in the lower troposphere below 5 km and in the tropopause region 9 km 
and above. These two turbulent regions were separated by a relatively quiescent gap of 
about 4 km. The strongest turbulence layers were well correlated with inflections in the 
temperature profile, thin stable or inverted temperature layers, and regions of high wind 
shear. A more detailed comparison of meteorological parameters with the C^(h) data 
appears in the next section. 
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Figure 3. Data from Thermosonde Flight 2, 5 May 1972. 

(a) Refractive-index structure coefficient (h) ascent data. (c) Low Richardson! number (Ri < 1,0). 

(b) Refractive-index structure coefficient C^(h) descent data. (d) Vector wind shear squared. 


(e) Temperature T and wind speed v versus altitude above MSL. 
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Figure 9. Data for Thermosonde Flight 9, 6 September 1972 
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APPUCATIQNS OF THERMOSONDE DATA 


Correlation with Meteorological Data 

Upper altitude wind and temperature information from the standard radiosonde flights were 
plotted in figures 3 to 9, in order to provide background meteorology and to indicate correla- 
tions with the microthermal data. The wind field is specified by flow speed and direction. 
Wind speed profiles appear directly in the figures. Wind direction data are used along with 
speed information to compute vertical wind shear profiles. The shear is computed by nu- 
merical differentiation of vector wind vertical profiles. The vertical spacing between each 
pair of wind data points in the radiosonde record is approximately 300 m. This is the resolu- 
tion element for shear computation. The magnitude of the vector difference between ad- 
jacent wind data points was divided by the 300-m interval to yield the vertical vector wind 
shear. A profile of the square of this shear appears in each figure. The radiosonde ambient 
temperature data were available only at relatively few significant levels, rather than the 
regular intervals of the wind data. Significant levels were chosen at those altitudes at which 
an inflection point occurred in the temperature vertical profile. Inflection points were stable, 
or inverted temperature lapse layers, or simple changes in temperature profile slope. 

Spikes or layers in the CjJ, (h) data appear to be well correlated with both temperature pro- 
file inflection points and vertical vector wind shear maxima. Good examples of this behavior 
are found at the 7.5- and 1 0.6-ktn levels in flight 2 and the 1 3.0-km level of flight 6. These 
correlations were achieved despite the separation of up to several hours and tens of kilo- 
meters between radiosonde and thermosonde flights. Often the temperature inflections and 
wind shear maxima were also well correlated with each other. 

Table 2 provides a semi-quantitative view of these correlations. All major features of the 
temperature profile were listed along with the altitude, or altitude interval of occurrence. 
Major features included stable or inverted temperature layers and inflection points in the 
lapse rate. The vector wind shear and temperature structure coefficient profile were then 
examined by the author and the presence or absence of major features (maxima) in these 
profiles were recorded in the tabic. The major features in all profiles w-ere rated as weak to 
moderate (o) or moderate to strong (•). In order for a major feature to be scored as 
present (correlated) with the temperature profile feature it had to be within ± 200 m of the 
temperature feature altitude in column 3 of table 2. Even the weak profile features were 
readily distinguishable by the height of their maxima above the average profile envelope. 

Note that CjJj(h) maxima are labeled with A or D for observations on ascent or descent, 
respectively. 
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Table 2 (Sheet 1) 

Correlation of Turbulence and Meteorological Events 


Flight 

No. 

Temperature 

Wind Shear 
Presence and 
Intensity 

Turbulence 

Profile 

Feature 

t 

Altitude or 
Altitude Range 
(km) 

Feature 

Intensity 

* 

Presence and 
Intensity 

Ascent (A) 
or 

Descent (D) 

2 

SL 

2.0- 2.4 

o 





SL 

3.8- 4.0 

m 

o 


A/D 


AT 

5.0 

0 





SL 

7.2- 7.6 

• 

# 

•/o 

A/D 


AT 

9.7 

o 

o 

• 

D 


IL 

10.6-10.7 

• 

• 

• 

D 


SL 

10.7-11.4 

• 

» 

• 

A 


SL 

11.8-12.9 

• 

• 

o 

D 


AT 

13.7 

• 




4 

AT 

1.5 

o 

o 




SL 

3.8- 4.6 

• 

o 

o 

D 


AT 

7.4 

o 





AT 

9.7 

o 

o 

0/0 

A/D 


IL 

11.5-12.2 

+ 

» 

o 

D 


AT 

13.9 

o 

o 




IL 

15.5-15.8 

• 

• 

o 

D 


IL 

16.3-16.6 

• 

• 

o 

A 


AT 

18.6 

o 

o 



5 

SL 

1.3- 1.6 

1 o 

o 

1 •/• 

A/D 


AT 

3.8 

' • i 


« 

A 


AT 

5.0 

° 


• , 

D 


AT 

6.3 

° 

o 

• 

D 


AT 

7.6 

o 





AT 

13.4 

0 

• 

• 

A 


IL 

13.9-14.4 

• 

• 

• 

A 


AT 

15.2 

• 


o 

D 


IL 

16.0-16.6 

# 





AT 

21.0 

o 




6 

IL 

2.2- 2.5 

« 

o 

• 

A 


AT 

3.2 

o 

o 




IL 

4.8- 5.1 

• 

o ^ 

• 

D 


f SL = stable temperature layer * o Weak to moderate 

IL = inverted temperature layer • Moderate to strong 

AT = change in temperature gradient 



Table 2 (Sheet 2) 

Correlation of Turbulence and Meteorological Events 
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Table 2 (Sheet 3) 

Correlation of Turbulence and Meteorological Events 


Flight 

No. 

Temperature 

Wind Shear 
Presence and 
Intensity 

... 

Turbulence 

Profile 

Feature 

t 

Altitude or 
Altitude Range 
(km) 

Feature 

Intensity 

* 

Presence and 
Intensity 

Ascent (A) 
or 

Descent (D) 


AT 

21.0 

o 

o 

o 

D 


IL 

22.2—23.5 

• 

o 

o 

A 

9 

AT 

1.2 

o 





SL 

1.2 - 3.1 

o 

o 


A/D 


AT 

5.1 

o 

o 




AT 

9.5 

o 

• 


A/D 


AT 

11.2 

# 

• 

•/o 

A/D 


AT 

12.0 

o 


o/o 

A/D 


AT 

13.4 

• 

# 




IL 

13.8-14.1 

• 

0 

o/o 

A/D 


IL 

17.0 

# 


o 

D 


AT 

18.6 

o 





AT 

24.2 

• 

o 




The degree of correlation among the columns of table 2 is quite good. About 75 percent of the 
temperature features are correlated with either or both wind shear and (h) maxima. The 
best one-to-one correlation, 73 percent, appears between (•) temperature features and wind 
shear maxima. The second best correlation of 7 1 percent is between (•) temperature features 
and C^(h) maxima. The presence of wind shear at the altitude of a (•) temperature feature is 
found to essentially guarantee (87 percent correlation) the appearance of a (o) or (•) turbu- 
lence maximum. Thus both temperature profile features and wind shear maxima are found to 
be good predictors of turbulence activity at high altitude. A closer look at the data of table 2 
reveals that the best correlations of the type noted above are found between 2 km and 10 km 
in the troposphere and at the tropopause. 

Wind shear and temperature information, when combined to form a vertical profile of gradient 
Richardson number, provide a truly quantitative estimate of atmospheric fluid stability, and 
thus can lead to a prediction of the presence of turbulence activity. The gradient Richardson 
number, 
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is defined in terms of g = 9.8 m s' 2 , temperature gradient 3T/3h, adiabatic lapse rate y a = 

+9.8 K. km' 1 , temperature T(K), and vertical vector wind shear 3V/3h . Regions where 
Ri < 1.0 are classed as unstable, while it is generally agreed that Ri <0.25 is required for tur- 
bulence initiation. Values of Ri <1.0 computed from the radiosonde data are indicated in 
figures 3 to 9 by the solid lines below the C^(h) data. Small values of Ri are achieved with 
small (3T/3h +7 a ), large i3\^/3h| 2 , or a combination of both. The concept of stable tempera- 
ture sheets (small negative 3T/3h) or inversion layers (positive 3T/3h) associated with tur- 
bulence seems contradictory to the attainment of low Ri numbers. Actually, as Scorer 1 2 has 
pointed out, initial tilting and thickening of the stable sheets could cause very local decreases 
in Ri numbers, which in turn would permit turbulence onset in an otherwise statically stable 
environment. Often the wind shear present is sufficient to produce a low Ri number for the 
sheet or structure as a whole. On this basis we would expect best correlation between small 
Ri and maxima of (h) in those regions where either data set is computed for atmospheric 
layers of relatively large thickness. This is because of the resolution of Cj^fh) data (75 m or 
1 50 m) and radiosonde data (300 m). Local decreases in Ri occurring in thin layers would 
not be resolved. Reiter and Lester 1 3 have demonstrated the striking dependence of calcu- 
lated Ri value on the altitude interval used in gradient calculations. From their evidence it is 
quite possible that if higher resolution meteorological data were available, the computed Ri 
could decrease to values less than 0.25 and predict directly the observed turbulence. As it 
is, the computed Ri values are typically large, rarely droping below 1 .0 or entering the critical 
region below 0.25. Regions where they do are well correlated with (h) maxima. Good 
examples of this behavior are found at the 1 1.1-km level of flight 5 and the 2.2-km level 
of flight 9. These could be examples of well developed turbulence where the Ri number is 
reduced for the structure as a whole. In other regions minima in Ri, while not less than 0.25, 
do occur at the same altitude as C 2 (h) peaks. These regions suggest the need for higher reso- 
lution meteorological data. 

The Cj^(h) data characteristics of strong thin layers, well correlated with the presence of 
temperature iiificetions, stable layers and wind shear, tend to support the explanation for 
high altitude turbulence generation in terms of breaking wave motion and Kelvin-Hclmholtz 
instability. The interested reader is referred to Scorer 1 2 and Woods 14 for readable accounts 
of the turbulence theory. The essential points of the theory are as follow's. Stable stratifica- 
tion or layers in the atmosphere occur at the boundary between air masses or as an upper 
limit for convective activity. Wind shear frequently is present or develops across this layer 
boundary. Internal gravity waves (IGW), induced by air flow over even minor terrain features, 
propagate continuously through the troposphere. Upon reaching higher altitudes and thus 
lower densities, the IGW increase in amplitude due to conservation of energy. Either the IGW 
or large vertical wind shear components tend to tilt the stable temperature layers in a periodic 
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fashion. Often the result of this tilt would be a simple sinusoidal shape for the stable layer. 

The tilted layer would typically have an initial thickness of a few tens of meters, amplitude 
and wavelength of hundreds of meters, and horizontal extent of several kilometers. Examples 
of this action appear in the periodic cloud patterns that occur at various levels in the tropo- 
sphere and are easily visible to the human eye at ground level. The clouds act as tracers to 
make visible the tilting phenomena which occurs predominately in clear air above convective 
activity. The roughly sinusoidal layers formed by tilting are called Kelvin-Helmholtz rolls 
or billows. While 1GW are always present in the atmosphere, they predominate in the lee of 
wind flow over mountains. Thus strong turbulence should be associated with mountainous 
terrain. Excellent support for turbulence in these mountain waves was found in the White 
Sands Missile Range thermosonde data reported by Bufton 8 . 

As the billows develop they tend to overturn and become dynamically unstable. The over- 
turning is analogous to the behavior of ocean waves approaching the shore. Wave troughs 
tend to recede opposite to the direction of mean flow, while wave crests are proceeding in 
the direction of mean flow. The result is a thickened layer where troughs and crests meet. 

As the wave proceeds the crests overturn, a curl develops, and the wave breaks into turbulent 
activity. It is the thickening of the wave as trough meets crest that produces instability. As 
Scorer 1 2 pointed out, R ; is inversly proportional to the square of layer thickness and wind 
shear. As a result very local decreases in R } are most likely in thick sections of the layer 
where shear is present, if R ; becomes less than 0.25, turbulence onset occurs. The turbulence 
expands to four or five times the original layer thickness, acquiring energy by feeding on the 
available wind shear. 

As a consequence of the turbulent mixing.two new stable layers are created, at the top and 
at the bottom of the original but now thickened and turbulent layer. As time progresses 
the two stable layers can also break down into turbulent activity, while the original turbulence 
dies away due to depletion of input wind shear. Good examples of this behavior appear to 
be present in the Cp(h) data of figures 6 and 9. In particular the two strong layers of C^fh) 
at 9.2 and 10.1 km in figure 9 on either side of the strong shear at 10.7 km may be a fine 
example of a later state of turbulence development. 

Comparisons with Stellar Observations 

The first flight series of the new thermosonde payload were accompanied by nearly simultaneous 
stellar observations with a ground-based telescope. The major motivation for this work was 
the desire to predict observed optical statistics by means of measured turbulence structure 
and current optical propagation theory. The side-by-side comparison of turbulence and 
optical data is illustrated in the experiment concept of figure 10. Although the optical 
propagation data were generated by observations of bright stars near zenith, the propagation 
link is identical to direct detection of down-link laser radiation from a source in space. The 
Cj^(h) profiles and comparisons with propagation statistics in this experiment serve to 
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STAR 



Figure 10. Experiment Concept for Comparison of Stellar and Turbulence Data 





illustrate the application of the thermosonde technique. Turbulence structure at the tropo- 
pause has traditionally been held to be the major cause of stellar scintillation. The evidence 
at hand both from radar studies and the present thermal sensor data suggest that while strong 
turbulence is often correlated with the tropopause, other regions, in particular, the lower 
troposphere, exhibit large values of Cj^(h)- The important question to resolve is whether the 
weight assigned to Cj^(h) versus h by the theory for prediction of scintillation statistics is 
necessary and sufficient to predict these same statistics. 

The stellar and turbulence data presented here were obtained on several clear nights during 
July, August, and September 1972 at the Goddard Optical Research Facility. The balloon 
flight plus simultaneous optical observations were conducted within 1 to 2 hours of the 
regularly scheduled National Weather Service radiosonde flight from Sterling, Virginia at 
00:15 GMT. The standard radiosonde flight provided temperature, pressure, and wind 
velocity profiles for interpretation of the turbulence and optical data. Each turbulence 
profile used in this analysis was a composite of descent data and ascent data. Ascent data 
were used only as necessary to Fill in the gap below 1 or 2 km caused by telemetry loss 
on descent. A preference for descent data results from the desire to eliminate wake effects, 
or false turbulence readings, from the data. 

Optical propagation data consisted of measurements of stellar scintillation and were obtained 
with a 50 centimeter (20-inch) reflector telescope with aperture stopped down to a clear 
17.7-centimeter aperture. An EMI 955 8B photomultiplier tube (PMT) was located in the 
telescope focal plane and provided the raw stellar data. An optical filter together with the 
PMT S-20 photocathode restricted the optical bandpass to 150 nm centered at 520 nm. 

Figure 1 1 is a block diagram of the scintillation monitor for both day and night operation. 

A more detailed description of the detector head is presented by Button and Genatt. 15 
Analog processing electronics were used to compute the fractional modulation of scintilla- 
tion (MOD), 

true rms irradiance 

MOD = ; (4) 

(average - background) irradiance 

The raw stellar data were also fed into a real time spectral analyzer and digital integrator. 
Spectral analysis was from 0.16 Hz to 200 Hz with a 1-Hz resolution. Typical data run length 
was two minutes. The results were the strength of irradiance fluctuations (scintillation) and 
their temporal spectra. 

Optical propagation theory for the stellar observation case is cast in terms of integrals over 
the C^j(h) turbulence structure. Thus the first step in comparison of the two data sets was 
computation of the required integrals. These integrals were computed numerically by sum- 
mation of C^(h) points separated by an altitude interval of 1 50 m. For the typical descent 
rate of 10 m s' 1 and data time constant of 5 seconds, each C* (h) point was an average over 
three vertical resolution elements. 

f i 
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Figure 11. Scintillation Monitor for Day and Night Observations 
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The first integral. 


I 


i 



c£C0 dh. 


( 5 ) 


is a linear weighting and sum over C^(h) at all altitudes from surface level, h Q , to the maxi- 
mum altitude, h, of reported data for the particular flight. An integral of this form appears 
in the theoretical description of optical wavefront statistics for plane wave propagation of 
starlight to a ground receiver. An estimate of Ij can be derived from the Tatarski 1 6 result 
for image motion variance, 

= Kj D“ 1/3 sec<?/ C 2 (h)dh, (6) 

with the additional assumptions that, 

8 = stellar zenith angle = 0° 

D = aperture diameter = 0.177 in 
K, < 2.91 

a 2 ^ = 0.5 arc sec 2 (typical result for 0 = 0° and D = 0. 1 77 m) 

The result is a prediction that 


I, > 4.7xl0" 13 m ,/3 

The I, integrals computed from the C 2 (h) data profiles are listed in table 3. The agreement 
between predicted and observed values is with a factor of 2 to 4. Contribution to I j from 
the surface boundary layer (first 10 m) ranged from 22 percent for flight 6, to 47 percent 
for flight 9. The factor of 2 to 4 agreement seems quite good considering the lack of direct 
measurement of and the uncertainty in the value of K ( from the theory. 

The irradiance statistic MOD is found in the wave optical treatment of stellar irradiance 
statistics of Tatarski.^ It can be squared and written as— 

MOD 2 = K 2 secT 0 I 2 , (7) 

for observations at stellar zenith angle 8 . Theory yields the value 7 - 3.0 for zenith angle 
dependence in the limit of a large telescope aperture. The experimental data of Protheroe 1 7 
can be used to estimate 7 = 2.5 for D = 0.1 77 m. I 2 is the second integral, 

/ ^max 

C 2 (h)h 5 / 6 dh, (8) 
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Table 3 


Comparison of Integrals Over C 2 (h) and Fractional Modulation Squared (MOD 2 ) 
Values as Computed from Thermal Sensor (Comp) and Optical Data (Data) 


Flight 

Number 

— 

h 

max 

(km) 

'*1 

(m ,/3 )X 10 12 

— 

h 

(m 7/6 )X 10 9 

, h 

(m 7/3 )X 10 5 

M °DL„ 

Zenith 

Angle 

(degrees) 

MOD 2 

Comp 

4 

18.3 

1.8 

1.1 

5.9 

0.(8 

37 

0.13 

6 

20.2 

1.0 

1.2 

63 

0.21 

34 

0.13 

7 

24.4 

1.4 

0.85 

3.8 

0.33 

66 

0.61 

9 

24.1 

1.8 

1.6 

1.7 



0.21 

14 

0.10 


over C 2 (h) turbulence structure weighted by the 5/6 power of altitude. The K 2 proportion- 
ality “constant” includes wavelength dependence, aperture averaging factor and numerical 
scaling. No attempt is made here to report or derive K 2 from the theory, since additional 
data on the jrradiance correlation length is required in order to accurately predict aperture 
averaging. Atmospheric dispersion effects discussed by Young 18 would also have to be 
taken into account. These effects influence the data because of the 1 50 nm optical bandpass 
and observation at angles other than zenith. Values for I 2 computed in a similar numerical 
fashion to those for are listed in table 3. Also included there are the corresponding 
M0D Data values and the zenith angle of observation. The average (1 .2 X 1 0 ~ 9 m' 7/6 ) of 
the four I 2 values compares roughly with a value of 6.0 X 10' 10 m _7/6 reported by Hulett 19 
on the basis of widely scattered scintillation measurements. The trend in 1 2 among the flights 
is repeated by the scintillation data with only little relative scatter. This argues for data 
validity. The departure of the two data sets from a one-to-one correspondence may be ex- 
plained by the aperture averaging and dispersion effects noted above. These effects are 
probably most severe for the flight 7 MOD 2 data taken at 9 = 66°. 

Data 

The temporal spectrum of the jrradiance pattern sampled by the telescope aperture provides 
another test of data and theory. Aperture averaging and zenith angle dependences appear 
clearly in the analytical result. In this test the C^(h) profiles plus the wind speed normal to 
the stellar line of sight are inserted into the theory and the result is compared with irradi- 
ance spectral density data. 

The normal wind speed component was obtained from the vector winds reported by radio- 
sonde. The conversion, as reported by Young, 18 


V^(h) - V Q (h)cos0 [1 + tan 2 0 sin 2 (fi - a)] 172 . 


(9) 
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employs wind speed V Q (h), wind direction a, stellar zenith angle 9, and stellar azimuth 

0 . 

The theory is derived according to the method of Tatarski. 16 The starting point is his equa- 
tion (13.38) for irradiance frequency spectrum in terms of the two-dimensional spatial 
spectrum of irradiance and a filter function for a ciruclar unobstructed aperture. The form 
below is obtained for the vertical propagation path: 



( 10 ) 


where: 

W, (f) = irradiance spectral density, 

V N — component of wind velocity normal to the propagation path, 
k = spatial frequency in aperture plane, 

f = temporal frequency, 

D = aperture diameter, and 

k Q = 2 ir[h 0 ~ optical propagation constant. 

This result states that the amplitude of irradiance spectral density at each frequency is a 
linear combination of (h) weighted by inverse wind speed and a sum over the spatial 
spectrum of turbulence, a diffraction term, and aperture filter function. 

The integrals of equation (10) were computed numerically on a digital computer after inser- 
tion of the observed C^(h) and V N (h) profiles, the 0.177 m aperture and the stellar zenith angle. 
An 8-point Gaussian quadrature formula was used to perform the inner integration. Typical 
integral limits for this procedure corresponded to a /c of 6.3 m -1 and of 57 m -1 for inner 
and outer limits, respectively. The limits h Q and h m aK appropriate for each flight were used 
for the outer integral. The height interval was again 1 50 m and the standard trapezoidal 
integration method was used. 

Irradiance spectra computed in this manner for flights 6 and 9 are indicated by the solid 
lines in figure 12. The open and closed circles indicate the corresponding stellar data. Each 
irradiance data spectrum was computed for a 2-minute segment of data used previously to 
compute the MOD£ Jta values reported in table 1 . Note that spectral density functions 
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normalized to unity are plotted. There is quite good agreement between data and prediction 
under different wind conditions of the two flights. The largest differences occur at low 
frequencies and for the lower wind speeds of flight 9. Possible error mechanisms at work 
here may be a breakdown in assumptions of equation (10), truncation of the Kolmogorov 
spectrum at low spatial frequencies, or too much weighting of high-altitude turbulence. 

The constants used to normalize the spectra of figure 1 2 were computed by integrating 
Wj(f) of equation (10) over the frequency range of interest. These constants are predictions 
of MOD 2 based only on the C^(h) and V N (h) profiles and theory, 

MODLw '/ f,l) 

*^0 

The upper limit f max was set at 200 Hz. Values of MOD 2 , omp are listed in table 3 for 
comparison with experimental MOD| ata . Agreement between predicted and experimental 
results is within a factor of two for all flights. The error between the results is a monotonic 
function of the zenith angle. This suggests the presence of an additional zenith angle depend- 
ence in the data, perhaps that due to the dispersion effect noted above. A more detailed 
discussion of these stellar comparisons with implications for the optical propagation theory 
is given by Button. 20 

Daytime Turbulence and Stellar Results 

Turbulence profiles have also been obtained during the daytime. Figure 8 is a C^(h) daytime 
profile. Note that while similar in structure for both day and night the turbulence layers 
appear stronger in the lower atmosphere during the day. This trend is substantiated by the 
work of other researchers. 1 > 2 Upper altitude turbulence also appears stronger during the day. 
Although the sharp inversion layer at 1 9 km suggests that strong turbulence should be found, 
it is not possible to be certain that these data are correct. Even if incorrect there would be 
little net effect on prediction of optical propagation statistics. Further evidence for the 
validity of these layers will have to wait for the availability of simultaneous optical data dur- 
ing the daytime. Some preliminary daytime stellar observations have been made with a 
modified version of the nighttime stellar scintillation monitor. They indicate a factor of two 
increase in scintillation strength (MOD) over average nighttime results. The signal-to-noise 
ratio for daytime stellar work was improved by a background discrimination technique out- 
lined in figure 11. The incoming stellar radiation and the sky background were chopped at 
a 2-kHz rate by clear and opaque sectors of a glass disk before detection by the PMT. The 
field of view was adjusted to just equal one chopping cycle on the glass disk. Thus the amount 
of background entering the PMT was constant in time despite the chopping action. The raw 
stellar data from the PMT were then bandpass-filtered to remove the dc background com- 
ponent. A standard AM demodulation technique of an absolute value circuit followed by a 
low-pass filter recovered the scintillation data. 
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The daytime thermosonde measurements have an additional temperature error source in the 
presence of direct solar radiation. The effect of solar radiation on the thermal sensor element 
itself is one possible source of AT errors. Solar heating of the probe support structure with 
transfer of heat to the proble is an additional factor. If both probes are heated equally, no 
AT error should result. A typical balloon flight, however, involves pendulum motions and 
payload spin. The results are continuous changes in the solar aspect angle and the possibility 
of occasional shading of one probe. From the calculations and data reported by Ney, Maas, 
and Huch 21 the solar radiation error for 2.5-pm diameter probe wire at high altitude can be 
estimated at several hundredths of a degree celsius. This can lead to appreciable error 
in the c£ (h) measurements. Radiation effects of this nature are more predominate at high 
altitude because of less efficient heat removal by convection. As a result the daytime indi- 
cation of strong turbulence near 20 km and above in figure 8 should be regarded with some 
suspicion. 

Prediction of Laser Propagation Data 

The major application of the microthermal turbulence data (C^(h)) to the prediction of laser 
propagation data is found in the downlink case. Laser radiation orginating in a satellite or 
spacecraft propagates downward through the entire atmosphere to a receiver station at ground 
level. The laser link would be used for high-data-rate communication. The ground receiver 
terminal would take one of two forms: a large photon bucket for direct detection of all 
possible incident laser radiation, or a heterodyne system for coherent detection. The space 
terminal is far enough removed that laser radiation is essentially plane wave by the time it 
reaches atmospheric turbulence. Thus the problem is directly analogous to the stellar ob- 
servations case described earlier. In direct detection the important information is the ir- 
radiance pattern sampled by the telescope aperture. Turbulence data enters through the 
theoretical prediction of irradiance statistics in terms of 

C 2 (h)h 5 ' 6 dh 

as in equation (8). In coherent detection, phase information is most important, and the 
propagation statistics contain integrals of the form 

C 2 (h) dh 

as in equation (5). A number of interesting propagation statistics can then be formed by 
inserting the above two integrals into the theory. The two most important statistics are 
aperture-averaged log-amplitude variance. 
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( 12 ) 


f oo 

c 2 = A(D,0,<?g 2 ) 0.56 k 7 / 6 sec# 11 ^ 16 / C 2 (h)h 5 / 6 dh, 

_/n 


where: 

A = aperture averaging factor, 

D = aperture diameter, 

6 = stellar zenith angle, 

k = optical wave number and 

h = altitude, 

for the direct detection case and coherence aperture diameter; and 


where: 

X 

e 

h 


r Q = 1.2 x 1 0~ 8 X 6 / 5 


[f 


C 2 (h)dh 


- 3/5 


(cos 0) 3/S , 


= wavelength in microns, 
= stellar zenith angle and 
= altitude, 


for the coherent case. 


(13) 


The uplink laser propagation case is also of interest for laser beacon and communication link 
possibilities. The primary statistic is received log-amplitude variance. Although the h 5/6 de- 
pendence is still present, the plane wave assumption is no longer valid. Account must be taken 
of transmitter beam size and divergence. Titterton 22 has indicated that the finite beam prop- 
erty may be used to reduce uplink scintillation by focusing the beam on the region of major 
turbulence in the upper atmosphere. Obviously the location and strength of turbulence 
layers are extremely important for this analysis. 

A further application involves resolution looking down through the atmosphere from a 
space platform. This has relevance to the areas of reconnaissance and remote sensing. Fried 23 
indicates that resolution should be proportional to the integral 



(14) 


This integral plus the two shown previously must be evaluated for quantitative prediction of 
the desired propagation statistics. In the present comparison of stellar observations and 
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turbulence data reported here, the many data points of each C^fh) profile were inserted and 
the integration performed numerically. It would be much simpler to assume a model for 
C^(h), perform the required integration analytically, and then adjust the model parameters 
to account for the time distribution and strength of turbulence. 

Numerous authors have suggested an exponential decrease in C^(h) with a separate maximum 
or bump to account for upper altitude effects. This in effect establishes an envelope for 
the behavior of multiple turbulence layers. As an example the profile could be written 


C$<h) - A 


exp (-h/H 0 ) 
h“ 


+ B exp 


(h-M 

20 2 


( 15 ) 


The parameters A, a, H 0 , B, £1 and Hj could be estimated from the thermosonde data by 
noting the relative strength of turbulence in adjacent atmospheric regions. A more tractable 
model might be a series of delta function (5) layers at regular height intervals (e.g., 1 km), 


C 2 (h) = Z VCh-hj). 
i 


(16) 


Delta function amplitude (A t ) would be adjusted to correspond to thermosonde data in the 
neighborhood of altitude h r The appropriate integrals are easily evaluated with delta func- 
tions in place of exponential and altitude power low dependencies. Fortunately altitude de- 
pendence of the propagation statistics is not critical on exact turbulence location. Accuracies 
within a km or so are sufficient. 


Application to CAT Detection 

In the current research C|(h) results are used to compute strength of refractive-index fluc- 
tuations C^(h) and thus the optical effects of atmospheric turbulence. These same tempera- 
ture fluctuations are, in fact, just evidence of density inhomogenities in the atmosphere. 
Since the atmosphere will not support local pressure gradients, warm parcels of air must be 
less dense than cold parcels at the same altitude. Thus temperature acts as a tracer for den- 
sity and C 2 (h) acts as a tracer for turbulence-induced density fluctuations. Density fluctua- 
tions in the range 20 m to 200 m form the clear air turbulence (CAT) that adversely effects 
aircraft operations. The small scale (of order 1 m) microthermal fluctuations measured by 
the thermosonde technique in terms of C 2 (h) depend for their existence on the cascade of 
turbulent energy from the larger scale structures. The correlation between small scale tur- 
bulence and large scale turbulence (CAT) has been demonstrated by recent radar studies of 
the upper atmosphere. In several cases 24 25 the radar observations (sensitive to turbulence 
dimensions of order 30 cm) have corresponded to simultaneous aircraft encounters with 
clear air turbulence. 
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In addition, the thermosonde data show consistent correlation between regions of high 
dynamic instability (Kelvin-Helmholtz billows) and C^.(h) peaks. The Kelvin-Helmholtz 
billows are considered a fundamental origin for CAT effects. The best application of the 
thermosonde technique to the CAT detection role would be, first of all, in side-by-side oper- 
ation with groundbased CAT radar systems. Then it may be desirable to launch thermo- 
sondes from standard radiosonde sites on a systematic basis to monitor a heavily traveled 
air corridor, and compare the data with aircraft turbulence encounters. Furthermore, once 
the spectral shape for turbulence wavenumbers is known, the strength of small scale turbu- 
lence can be used to predict quantitatively the strength of CAT. There is considerable exper- 
imental evidence to suggest that the predominate spectral shape is the (-5/3) power of wave- 
number predicted by the Obukhov-Kolmogorov theory for the inertial subrange. 


COMMENTS 


The micro thermal profile measurements reported here suggest that upper altitude turbulence 
structure is concentrated in various layers. While the strongest layers of turbulence were 
found at the surface boundary layer and associated with temperature inflections and shear 
layers in lower troposphere, significant contributions to C 7 (h) were recorded from the 
tropopause region and occasionally from the stratosphere. It is the strong correlations with 
meterological parameters and the predictable cascade of energy from large to small structures 
that suggest the application of the thermosonde to CAT detection. The observed turbulence 
combined with current theory works well to predict not only strength of scintillation, but also 
the observed temporal spectra, if wind velocity profile data are included. Agreement is good 
between prediction and data despite the expected non-stationarity of turbulence. The present 
method of differential thermal sensor measurements from a balloon payload yields vertical 
cross sections of a C 2 (h) structure that varies widely in time and space. Because of the hori- 
zontal sheet nature of this structure the thermal sensor data should be a reasonable first 
estimate of Cj^(h) strength and distribution. If the application is limited to comparison 
with simultaneous optical propagation data, results. within a factor of two should be expected. 
More accurate comparisons, prediction, and tests of propagation theory will require three- 
dimensional measures of turbulence in addition to the profile samples. 

The most pressing need in future investigations is a measurement program for turbulence 
spatial spectrum in the upper atmosphere. Until the presence of an inertial subrange and the 
validity of Obukhov-Kolmogorov statistics is verified, few quantitative comparisons can be 
accepted. It will also be necessary to define the limits of non-stationarity of upper altitude 
turbulence in order to extrapolate the point thermal measurements. Finally the same com- 
parisons of turbulence and optical data should be made during daylight hours to assess the 
effect of diurnal variations in turbulence strength at high altitude. 


Since the first publication of this material in November, 1973, as a GSFC X-document2 6 
two extremely important and related turbulence measurement programs have been reported 
in the literature. These independent programs are the work of Coulman 27 and Barletti 
et al. 2 8 The reported measurements are quite similar to those of the thermosonde device, 
in that micr o therma 1-sen sor equipped radiosondes were used to obtain the turbulence profile 
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data. The major significance of these efforts is the direct telemetry of raw turbulence data 
to the ground for subsequent data analysis. The latter two references should be consulted 
for interpretation of and additional insight into the results reported here. 
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